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I. INTRODUCTION 


The theory of open quantum systems is well estab¬ 
lished when a Markovian approximation applies [l|. In 
this situation the density matrix evolution is defined by 
a Kossakovski-Lindblad equation. In addition, the sys¬ 
tem dynamics can be read in terms of an ensemble of 
stochastic trajectories developing in the system Hilbert 
space [H-Q- Different stochastic Schrodinger equations 
characterize the ensemble dynamics (unravelings). The 
diffusive case 0 corresponds to multiplicative Gaussian 
white noises. It allows to describe quantum systems sub¬ 
jected to a continuous measurement process [2-13 
as to formulate dynamical wave vector collapse models 


In the last years, with the goal of establishing a non- 
Markovian extension of the standard Markovian open 
quantum system theory, different research lines were 
opened. In particular, the stochastic ensemble represen¬ 
tation of a quantum system coupled to a bosonic bath 
0 ini triggered the study of stochastic Schrodinger 
equations driven by multiplicative non-white Gaussian 
noises. As demonstrated in the seminal contributions 
of Strunz and Diosi 0, the emerging time-evolution- 
equation involves a functional derivative of the wave 
vector. Due to this feature, the formulation of differ¬ 
ent derivations and perturbation schemes, which lead to 
evolutions without involving an undetermined functional 
derivative, is a topic that is of interest up to present time 
[l^ - [l8j| - On the other hand, exact closed expressions for 
the functional derivative were found in different physical 
situations (l9l - [^ . 

The non-Markovian quantum diffusion model [isini 
provided an alternative context for discuss]^ non- 
Markovian continuous measurement theory [25l - l27l] , non- 
Markovian extensions of spontaneous wave function col¬ 


lapse models [H, [l^] , and operator correlations beyond 
the quantum regression theorem . The formalism was 
also used to characterize specific physica l sy stems such as 
for example quantum Brownian motion [3ll j , complex op¬ 
tical arranges [s^, 113 , charge transport in organic crys¬ 
tals [131 and many body systems [ssj - 

A relevant advance in the field was introduced recently 
by Diosi and Ferialdi [13 • On the basis of a path inte¬ 
gral approach, they defined a generalized class of non- 
Markovian Gaussian stochastic Schrodinger equations. 
Similarly to the Markovian case, the generalized unravel¬ 
ing is parametrized by a set of complex noise correlations 
that do not affect the density matrix evolution [3,13 • For 
Markovian unravelings the degree of freedom introduced 
by the extra correlations can be set from symmetries con¬ 
straints such as the invariance of the unraveling under lin¬ 
ear system operator transformations under which the cor¬ 
responding master equation is invariant [^-Q- Fhe main 
goal of this paper is to develop similar symmetry analysis 
for the generalized non-Markoyian unrayeling [36l j and to 
find which constraints on the noise correlations arise. 

We show that both the generalized non-Markoyian 
Schrodinger equation and its associated density matrix 
eyolution always share the same symmetry property, that 
is, they are inyariant under arbitrary unitary changes of 
the system operator base. Therefore, in contrast to the 
Markovian case, we focus on which constraints on the 
noise correlations may arise when mapping the unrav¬ 
elling invariance with that of a microscopic system-bath 
Hamiltonian description. Different solutions to this prob¬ 
lem are found, which are based on different kind of map¬ 
ping between the noise and bath operator correlations. 
The emerging conditions under which the mapping is 
consistent are studied for a single dissipative channel dy¬ 
namics. Furthermore, the stretched relation of the gen¬ 
eralized non-Markovian unraveling with previous results 
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derived from quantum-measurement theory are ana¬ 
lyzed from this perspective. 

We base the present analysis on a generalized non- 
Mar kovian Gaussian Schrodinger equation that, in con¬ 
trast to Ref. [ 3 ^ . become written in an arbitrary non- 
Hermitian base of system coupling operators. Instead 
of a path integral formalism, here the wave vector evo¬ 
lution [Eq. ([T^ ] is derived by postulating a stochastic 
density matrix dynamics driven by multiplicative non¬ 
white Gaussian noises, where an undetermined contribu¬ 
tion is obtained from a condition of average trace conser¬ 
vation ( 1 ^ . 

The paper is organized as follows. In Sec. II, based on 
the stochastic approach outlined previously, we obtain 
the general non-Markovian Gaussian diffusive unravel¬ 
ing. The main results are obtained in Sec. Ill, where 
its symmetry properties are analyzed. In Sec. IV we ex¬ 
emplify the main results by studying the case of a single 
Hermitian coupling channel. Previous results obtained 
from quantum measurement theory are also studied and 
derived by using the present approach. The Conclusions 
are given in Sec. V. In the Appendix we define the main 
properties of the noise correlations as well as a general¬ 
ization of Novikov theorem valid for arbitrary mul¬ 
tiplicative complex Gaussian noises. 


II. GENERAL STOCHASTIC GAUSSIAN 
DYNAMICS 

Both in the Markovian and non-Markovian regimes, 
stochastic Schrodinger equations define an ensemble of 
system states, which can be written in terms of a stochas¬ 
tic density matrix Pst{t)- Its average over the ensemble 
of realizations, denoted as (•••), gives the system state, 
p{t) = {pstit)). In the present approach, we start by pos¬ 
tulating the time evolution of Pst(t), 

= -i[Hs,Pst{t)]-l{[pst] 

( 1 ) 

Here, Hs is the system Hamiltonian. A is coupling pa¬ 
rameter. The fluctuation operator reads 

T{t)=J2^o{t)L^, (2) 

a 

where {La} is the set of coupling system operators. The 
set of (multiplicative) complex noises {Za{t)} are Gaus¬ 
sian with null mean value {Za {t)) =0, and correlations 

Xap{t,s) = {Za{t)zp{s)), (3) 

and also 

■nafi{t,s) = {Za{t)zp{s)). (4) 

Notice that in contrast with previous analysis M and 
similarly to Ref. [3^ . here we are taking into account 


the correlations rjapit, s). These objects allow to consider 
arbitrary correlated noises (see Appendix). While the fol¬ 
lowing derivation is similar to that when rjap{t,s) —>■ 0 
here we show it in order to enlighten the main as¬ 
sumptions over which it relies. 

The unknown functional U[pst\ is determined by aver¬ 
age trace condition 

^Tr(pst(t)) = 0, (5) 

and the separability condition (pure state unraveling) 

PstW = IV'(i)) (V'(i)l- (6) 


This last requirement allow us to defining a stochas¬ 
tic Schrodinger equation for the system state IV' (t)) • In 
addition, this requisite implies that the stochastic map 
Pst{0) — t Pst{t) is a completely positive one [l|. There¬ 
fore, Eq. ([5]) is a sufficient condition that guarantees the 
completely positive property of the density matrix evo¬ 
lution p{0) —!> p{t). 

By averaging Eq. m and imposing condition ([S]) we 
get 

Tr {U[p,t]) = -iATr {{T (t) pst{t)) - {pst{t)T'^ (t))) . 

( 7 ) 

Novikov theorem gives an exact (functional) expres¬ 
sion for the average of a product between a Gaussian 
noise and an arbitrary functional of it. In the Appendix 
we provide its generalization for arbitrary complex noises. 
Hence, from Eqs. ([2]) and (IA.3I) the previous contribu¬ 
tions can be written as 


{t) pst{t)) = dsXap{t,s)La(^j^^'^ 

+ J dsr]ap{t,s)La(^^^^^y ( 8 ) 

and similarly 


The convention of sum over repeated subindices applies 
whenever the summatory symbol is not written. Given 
the commutativity of the trace operation, the indetermi¬ 
nation given by Eq. © can be surpassed by demanding 
the separability condition ([6]), which due the analytical 
property of the wave vector |V'(i)) implies 


(V^(0l 

5z*p{s) ■ 

. . 

Therefore, the contribution lA[pst] is determined in an 
unique way, 


dpstjt) ^ d\'tl){t)) 
dz/ 3 (s) 5zp{s) 




Spstjt) 

SzUs) 


= IV'W) 


W[pst]=zA / ds{Llxap{t,s)-Lar]al3{t,s)) 
Jo 


Spstjt) 

5zp{s) 


+H.C., 

( 11 ) 
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which in turn lead to the general non-Markovian Gaus¬ 
sian stochastic Schrodinger equation 

A = _iHs \'iP{t)) - i\F (<) IV'(t)) (12) 

-iX J ds{Lixa/3{t,s) - Lar]a/3{t,s))j^^^. 

This is the main result of this section. We notice that 
under the replacement —>■ L^, this equation recovers 
the stochastic dynamics of Ref. [s^. As demonstrated 
in the next section, this replacement is equivalent to a 
unitary change of the base of system operators. Hence, 
Eq. (1121) and the evolution obtained in |36l| are equivalent 
dynamics expressed in different operator bases. On the 
other hand, when rjapit^s) —^ 0 it follows the standard 
non-Markovian quantum diffusion model EH. In the 
next section, we also show that the underlying symme¬ 
tries of the generalized time evolution can be inferred 
and enlighten from Eq. (fT^ (Sec. III). In addition, this 
equivalent evolution allow us to recover in a simple way 
previous generalized stochastic dynamics obtained from 
quantum measurement theory |25l | (Sec. IV). 

Functional structure and density matrix evolution 


result that follows from Eq. (HU. 

We notice that Eq. dn does not depend explicitly on 
the correlations rjap{t,s). By performing a recursive ex¬ 
pansion in the interaction strength parameter A 0 , it is 
also possible to demonstrate that the density matrix evo¬ 
lution, order by order in A, does not depend of riap(t, s). 
Therefore, consistently with Ref. 0, we conclude that 
these correlations only modify the wave vector evolution. 
Their physical role is investigated through the following 
symmetry analysis. 


III. SYMMETRIES 

Here, we analyze the symmetries of the generalized 
stochastic Schrodinger evolution defined by Eq. ©■ 


A. Hermitian fluctuations 

When the fluctuation operator ([2]) is Hermitian 

(16) 

it follows the equality 


As in the standard case 00, we notice that Eq. 
m depends on the functional derivative of the wave 
vector. Most of the achievements performed in the last 
years 0 - 0 | can be applied to the generalized unrav¬ 
eling. In particular, by following the calculation steps 
performed in Ref. [ij| it follows 

=-iXGstit,s)Li3\'ijj{s)), (13) 


where the (functional) propagator is G St {t, s) — 


exp -i J* drTst (t) 


. Here, [• • •] denotes a time order¬ 


ing operation while the functional generator Tst (t) de¬ 
fines the time evolution m, that is, (d/dt) \ip{t)) = 
Tst (t) \tp{t)) . As demonstrated in Ref. 0 these ex¬ 
pressions allow us to perform consistent approximations 
(in the interaction parameter or in the noise correlation 
times) of the stochastic wave vector in both a time con¬ 
voluted and time convolutionless schemes. 

By introducing Eq. m in the stochastic evolution 
(HD and after performing the average over realizations [see 
Eqs. ([5]) and ([H])], we obtain the density matrix evolution 


c/ ^ 

■^p{t) = -i[Hs,p{t)]-iX J ds|xa/ 3 (t,s) (14) 



where the functional derivative is given by 

= -'i>'Gst {t, s) Lppst {s)GI^ {t, s ), (15) 

oz^(s) 




This relation follows by multiplying Eq. ([16]) by Z/ 3 (s) 
and using the definition of the noise correlations, Eqs. o 
and O- We notice that the weight of the functional con¬ 
tribution in Eq. m is precisely the difference between 
the two terms in the equality (ED- Hence, it consistently 
vanishes for Hermitian fluctuations. 


B. Invariance under unitary changes of system 
operator base 

In the derivation of Sec. H the fluctuation operator 
T (t) was written in the base defined by the set of op¬ 
erators {La}. Here, we explore which structure assume 
the evolution (I12|) when introducing a new base of op¬ 
erators {A^} related to the previous one by a unitary 
transformation. Hence, we write 

Af_i ='^^LaUafi, Lg = Afj^Uj^g, (18) 

a. ^ 

where the unitary operator U satisfies = Sap 

and This change of base allows to 

define a new set of noises {r^(t)} related to the original 
one as 

*"At(^) = ^ ^ Za(t) = ^ ( Ug^T^it), (19) 

Q! fl 
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such that the fluctuation operator can be rewritten as 

= ( 20 ) 

OL ^ 

The variational derivative in Eq. dm can be written as 

<^(') _E]_ EEl — r/t ‘^(') 

Szpis) ^Jo 6r^{T) Szfi{s) ^ ''^6r^{s)' 

( 21 ) 

where we have used that 8 r^{T)/5zp{s) = S(t — s)U^^ 
[Eq. (ITOl) ]. Hence, from the previous expressions, in the 
new base the general stochastic Schrodinger equation be¬ 
comes 

I m)) = -^Hs \m) - *A {t) \m) m) 

where we have defined the correlations 
Xfj.i'it, s) = E ^«MXa/3(<, s)uli^ = (r* (t) (s)), (23) 

a/3 

and also 

s) = ^ s)Ulp = {r^ (t) (s)). (24) 

aj3 

Equation (1221) and the previous two expressions show 
that the general evolution dHD is invariant under the 
unitary changes defined by Eqs. m and (HU). One can 
always choose an Hermitian base of operators {H.^} = 
(Hj)}, which explicitly demonstrate that Eq. (IT^ and 
Eq. (27) in Ref. [s^ are related by a unitary change of 
operator bases. On the other hand, from Eq. m it is 
possible to demonstrate that the density matrix evolution 
da is also invariant under the transformation (jlSp and 
da, that is, it can be rewritten in terms of the operators 
Afj^, the derivative 5pst{t)/5r^ {s), and the correlations 
Xtiu{t,s) [Eq. da]- 

Given that x^j.v{t,s) and r]^^(t,s) only depend on the 
correlations Xap{i,s) and r]ai 3 {t,s) respectively, the in¬ 
dependence of p{t) with respect to pa 0 {t,s) [Eq. (IT4l) ] 
also implies its independence with respect to r]^,v(t, s) in 
the new base. Therefore, this property is valid in any 
operator basis, conclusion consistent with the results of 
Ref. [13. 

White noises 

For white correlated noises 

Xap{t,s) = 8{t-s)5ap, T]apit,s) = 6{t-s)Cai3, (25) 

where Cap are complex coefficients, by using Eqs. (fTdll 
and da, the density matrix evolution becomes a Lind- 
blad equation 

^ = -i[Hs,p{t)] + X^i[L^,p{t)Li] + [LMt),Li]). 


Under the unitary transformation da, this equation re¬ 
mains invariant, that is, the only change corresponds to 
the replacements —>■ On the other hand, in ad¬ 

dition to this change, Eq. da results defined by the 
correlations [see Eqs. (l23l) and (l24l) ] 

= 5{t-s)5^v, (26a) 

Thi<^{t,s) = 8{t-s)^Ul^CapUlp. (26b) 

a/3 

While the master equation is invariant (does not de¬ 
pend explicitly) under the unitary transformation U, the 
stochastic evolution depends on it through the correla¬ 
tion Therefore, if one demand that both the 

master equation and the stochastic Schrodinger equation 
must have the same dependence on U (symmetry) it fol¬ 
lows that Cap —S' 0. Hence, r]ap{t,s) —)> 0. Based on a 
different approach, this result was developed in Ref. Q- 
Now, we ask if this kind of arguments are also applicable 
for the generalized non-Markovian Schrodinger equation. 

C. Mapping with microscopic Hamiltonian 
symmetries 

The previous analysis (valid for white noises) does not 
apply in the present framework. In fact, a redefinition 
of the noises allowed us to conclude that the general¬ 
ized Schrodinger dynamics dED and the master equation 
(1141) share the same invariance symmetry property un¬ 
der unitary changes of the operator bases. Thus, any 
constraint applicable to the noise correlations has to be 
based on a more fundamental requirement. With this 
motivation, we ask if the stochastic Schrodinger equa¬ 
tion has the same invariance symmetry property than a 
microscopic bosonic dynamics able to induce the same 
system dynamics. 

The total microscopic system-reservoir Hamiltonian 
Ht reads 

= + + (27) 

a 

where as before Hs is the system Hamiltonian, Hb is 
the bath Hamiltonian and the remaining contribution 
gives their interaction. Similarly, we introduce a unitary 
change of the system operator base [Eq. (IT^ ] 

= E L<.Ua^, La=J2 (28) 

a ^ 

and a new set of bath operators 

=''^Ua^R^i, (29) 

a ^ 

in such a way that the interaction contribution Hj can 
be written as 

Hi = ^La® Za = ^A^®R^,. (30) 

a fi 
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When tracing out the bath dynamics, p[t) = 
TrB[exp(—exp(+iti77’)], where TrB[- • •] is 
a trace operation in the bath Hilbert space and ps is the 
environment stationary state, given that TtbIpeZoi (t)] = 
TrB[pB exp(+ztiJB)Za exp(—itiJs)] = 0, the relevant 
statistical objects are the bath operator correlations (in 
an interaction picture with respect to Hb) n There¬ 
fore, we search under which conditions both noise corre¬ 
lations Xa /3 [t, s) and rjapit, s) can be related or mapped 
with the quantum bath correlations. 

Some conditions are imposed over the bath-noise cor¬ 
relation mapping. 

(i) The map has to be invariant under unitary changes 
of the system operator base. 

(ii) The density matrix evolution obtained from the 
stochastic and microscopic Hamiltonian approaches must 
be the same. 

(iii) In addition, the mapping must be consistent, that 
is, the resulting noise correlation matrix has to be posi¬ 
tive defined [see Eq. (IA.4I) in the Appendix]. 

Under the previous conditions, different mapping with 
different physical motivations can be proposed. A di¬ 
agonal map arises from a direct mapping between the 
stochastic and microscopic dynamics, which lead to the 
condition r]ai 3 {t,s) —>■ 0. A non-diagonal map is moti¬ 
vated from quantum measurement theory, which allows 
us to raise the previous condition, ijapitjs) yf 0, if some 
constraints are fulfilled. The diagonal mapping can be 
read as a particular case of the non-diagonal one. Nev¬ 
ertheless, for clarity, below each case is presented in a 
separate way. 

1. Diagonal correlation mapping 

The stochastic and microscopic density matrix dynam¬ 
ics can be put in one-to-one correspondence [see Eqs. o 
and dsni] under the following correlation associations 

Xap(t,s) ^ TTB[pBZlit)Zp{s)], (31a) 

ria.l3{t,s) O- TlBlPBZait) Zp{s)]. (31b) 

This correlation map is invariant under unitary changes 

of system operator base [condition (i)]. In fact, from Eqs. 
(HU) and dMl), in the new base Eq. m becomes 

XMi^(i,s) TrB[pBRl{t) R„is)], (32a) 

-H- TTB[pBRfi{t)Ry{s)\, (32b) 

which corresponds to the same mapping in the new noise 
and bath operator bases. 

Given that the average density matrix does not depend 
on the correlations ijapitjs), the equality Xapi^^s) = 
TrB[/9BZf (t) Zp{s)] guarantees condition (ii), that is, the 
stochastic and microscopic approach lead to the same 
density matrix evolution (see Ref. H). Notice that due 
to the Gaussian and bosonic properties of the fluctuation 
operator and reservoir respectively, it is not necessary to 
define the map for higher correlations. 


The map (1311) should allows us to associate a micro¬ 
scopic origin to both noise correlations Xap (t, s) and 
ilap{t,s). Nevertheless, by taking the relations (I5T1) as 
valid equalities, by using that i?j(t) = from Eq. 

(HB one immediately deduce (in an interaction repre¬ 
sentation) that the fluctuation operator T{t) should be 
Hermitian, implying the cancellation of the variational 
contribution in the stochastic Schrodinger evolution m- 
Hermitian fluctuations can only be mapped with a quan¬ 
tum reservoir at infinite temperature |14l [ 3 ^ . Given an 
environment at finite temperature, this contradiction can 
only be surpassed by imposing the correlation map 

Xapit^s) = TTB[pBZl{t) Zp{s)], (33) 

and demanding 

r]ai3it,s) =0. (34) 

In the Appendix we demonstrate that, without imposing 
any condition on the bath properties, these correlations 
satisfy condition (iii). Therefore, Eq. ((33l) and (IMt define 
a consistent correlation map which in turn recovers the 
standard non-Markovian quantum diffusion model [T^ . 
ini (Eq. (IT^ with rjapitjs) 0). This is one of the 
main results of this section. Gonsistently, Eq. dSl also 
implies s) = 0. 

The previous correlation map [Eqs. (l33l) and (l34l) ] 
guaranties that the stochastic and microscopic evolu¬ 
tion (at any temperature) have the same invariance 
symmetry property under arbitrary unitary changes 
of the system operator base. Notice that the alter¬ 
native mapping 77^/3 (t, s) = ^^[pbZq, (t) Zp^s)] with 
Xa/3 (tj s) = 0 is not consistent with condition (iii) [see 
Eq. (IA.4I) ]. In fact, for arbitrary complex noises it fol¬ 
lows Xa/3 (t, s) \t=s,a=P > 0. 


2. Non-diagonal (measurement-like) correlation mapping 

In diverse quantum optical arranges, the Markovian 
dynamics of a system is inferred from a measurement 
process performed on the environment. Different mea¬ 
surement schemes are defined by the quadratures of the 
bath, which in turn define different stochastic unravelings 
ii with the same average dynamics. In the present 
context, a similar degree of freedom can be introduced 
by defining a set of quadrature-like ii bath operators 
{Za{t)}, which read 

Z„(t) =^M„,,Z,,(t). (35) 

a' 

Here Maa' is an arbitrary complex matrix. When 
Maa' = Saa' the interaction bath operators [see Eq. 
are recovered. Therefore, consistently with the previous 
diagonal case, maintaining the base of system operators 
{La}, instead of Eq. (I5T1) . we introduce the generalized 
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correlation mapping 

Xai3{t,s) = TrB[pB2l{t)Z/3{s)], (36a) 

r]a/3{t,s) = TrB[pBZa{t) Zi3{s)]. (36b) 

The non-diagonal associations p6p are also invariant 
under changes of the system operator base [condition (i)]. 
In fact, from the Schrodinger evolution (1^ . Eq. (^5)) and 
(EH), it follows Eq. dSH) with the replacement (t) —)■ 
TZ^ (t) , where 

(t) = (t) = (t) . (37) 

Hence, invariance is guaranteed by a unitary transforma¬ 
tion of the matrix Maa' ■ 

Condition (ii) implies that the density matrix evolu¬ 
tion has to remains the same than in the diagonal case, 
which in fact corresponds to the dynamics derived from 
the microscopic Hamiltonian dynamics. Given that in 
Eqs. (IHTI) and (1551) the base of system operators is the 
same, and given that the density matrix evolution da 
does not depend on the correlations Papit) [l^, it follows 
the condition 

TtbIpbZI (t) Zp(s)] = TtbIpbZI (t) Zp{s)]. (38) 

Therefore, the bath operators {Za(t)} and the 
quadrature-like bath operators {Za (t)} must have the 
same correlations. Erom Eq. (1351) this condition explic¬ 
itly reads 

TtbIpbZI (t) Zp{s)] = Ml^Mpp,TrB[pBZl (t) Zp,{s)]. 

Even when the previous constraint is fulfilled, it is not 
possible to guaranty that condition (iii) is satisfied. That 
is, it is not possible to know in general if the quantum 
correlations defined by {Za (t)} [Eq. (1551) ] satisfy or not 
the positivity condition given by Eq. (1A.41) . 

In general, it may be difficult to satisfy the previ¬ 
ous requirements without imposing some condition on 
the matrix Map or alternatively on the properties of 
the bath. Notice that in the white noise approxima¬ 
tion TvB[pBZ'l{t) Zp{s)\ = 5{t - s)5ap [Eq. (ES])], 
Eq. (1551) is automatically satisfied for unitary matrixes, 
Ml^i^Mpa' = Sap (Stratonovich calculus). 

Assuming the consistence of the map (1361) with condi¬ 
tions (ii) and (iii), in contrast to the diagonal case, here 
the Hermiticity of the interaction Hamiltonian Hj (t) and 
Eq. (ITTl) do not lead in general to any contradiction. 
Therefore, given that Eqs. dMD are consistent with the 
correlations of a set of complex noises, and given that Eq. 
(1551) is satisfied, the non diagonal mapping [Eq. (1551) ] al¬ 
low us to define the correlations of the generalized Eq. 
EU from the microscopic underlying evolution. This is 
the second main result of this section. In fact, we ob¬ 
tained conditions under which the different unravelings 
related to Eq. © can be associated to the microscopic 
Hamiltonian dynamics. At this stage, notice that the di¬ 
agonal case can be read as a particular case of the present 
one. 


IV. EXAMPLES 

In the next examples we analyze the consequences 
of the previous constraints when applied to different 
system-bath interaction structures. The relation of the 
generalized stochastic unraveling with previous results 
based on quantum-measurement theory is also re¬ 
visited. 


A. Hermitian single channel 

The general evolution m takes a simpler form when 
one take into account only one single channel, T (t) = 
z{t)L, defined by an Hermitian operator, = L. There¬ 
fore, it follows 

= -iHs\tp{t)) - iXz{t)L\ilj{t)) (39) 

-iXL f 

Jo Sz[s) 

From a phenomenological point of view, one can propose 
different noise correlations x(t, s) and p{t^ s), as for exam¬ 
ple exponential ones. In this case, it is possible to explic¬ 
itly show [see Appendix, Eq. (IA.9|) ] that the extra cor¬ 
relation p{t, s) allows to smoothly change the stochastic 
unraveling between the standard non-Markovian quan¬ 
tum diffusion model [pit, s) = 0] and a pure stochastic 
Hamiltonian [pit,s) = x(t, s)]. On the other hand, the 
consistency of this evolution with the invariance symme¬ 
try of a microscopic Hamiltonian description is analyzed 
below. 

Eq. (I59t should be associated with an interaction 
Hamiltonian of the form H[ = L(g> Z, where the bath op¬ 
erator is Hermitian, Z = Z\ with (complex) correlation 
Mb[pbZ (t) Zis)]. If one try to impose the associations 
(1511) . it follows that xit^s) = p{t,s). This equality can 
only be satisfied by a real noise z{t) with a real correla¬ 
tion [see Eq. (IA.4I) ] . Consequently (t) = R (t) . Given 
that xitjs) = X*(^i'S) = vit^s) = p*it,s), it also implies 
Mb[pbZ (t) Z{s)] = TtbIpbZ {s) Zit)]. In general, this 
equality can only be satisfied when ps is proportional to 
the identity matrix, that is, a thermal state at infinite 
temperature. This contradiction, which also implies the 
cancellation of the functional contribution in Eq. (I39L 
is raised up by taking x(^! s) = Mb[pbZ (t) Zis)] and 
pit, s) = 0 , that is, consitently we recover the diagonal 
correlation mapping defined by Eqs. © and El re¬ 
spectively. 

Given that the interaction Hamiltonian only involves 
one single Hermitian operator Z, here it is not possible 
to introduce a non-diagonal correlation mapping El- 
Therefore, the correlation pit, s) in Eq. (1551) cannot be 
related to the microscopic interaction if the same invari¬ 
ance symmetry property is demanded. 
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Dephasing dynamics 

Dephasing dynamics is one particular physical case of 
the previous situation. It emerges when the (Hermitian) 
coupling operator commutates with the system Hamilto¬ 
nian, [Hs,L] — 0. Hence, in Eq. (IMl) . the variational 
derivative can be written as 

= -i\L m)) , (40) 


Noticing that {Za{t)} = {Z^{t), Z{t)} characterize the 
interaction Hamiltonian in the base of system operators 
{La} = {iL, —iL^}, demanding the positivity of the noise 
correlation matrix [see Eq. (IA.4I) ]. it is simple to check 
that the diagonal correlation mapping defined by Eqs. 
([33|) and (|34|) can only be consistently satisfied by intro¬ 
ducing one single complex noise Zc(s), T (t) = Zc{t){iL), 
such that Xc{t,s) = {z*{t)zcis)) = TTB[pBZ{t)Z'< (s)] = 
Z/c and r]c{t,s) = {zc{t)Zc{s)) = 0. 

Hence, Eq. (IT^ becomes [H{t) = 


result that follows straightforward from Eq. (fT51) . On 
the other hand, from Eqs. m and m it follows the 
density matrix evolution 

^Pit) = -i[Hs,p{t)]-X^ J dsK{t,s)[L,[L,p{t)]], (41) 

where K{t,s) = [x(^j s) -|-x*(t, s)]/2. We notice that this 
time evolution also arises from a pure stochastic Hamil¬ 
tonian (see Eq. (40) in Ref. [^) defined with a real 
Gaussian noise with correlation K(t, s). Therefore, all the 
set of ensembles defined by Eq. (l39ll and (IdHI) are equiv¬ 
alent in average to a stochastic Hamiltonian. This prop¬ 
erty is stretchy related to the dephasing property of the 
underlying fluctuations. 


B. Quantum optical-like microscopic interaction 

Eor quantum optical systems, where a Markovian ap¬ 
proximation applies, the correlations rjapit, s) (delta cor¬ 
related) define different ensemble dynamics associated to 
different measurement processes performed over the en¬ 
vironmental degrees of freedom [8|. Based on quantum 
measurement theory, Gambetta and Wiseman demon¬ 
strated that a similar result is valid in the non-Markovian 
regime Below, we show that those results can be 

recovered as diagonal and non-diagonal correlation map¬ 
pings associated to the same microscopic Hamiltonian. 
In fact, the degree of freedom associated to the different 
measurement processes here can be related to alternative 
definitions of the bath operators {Za{t)} [Eq. (1551) ]. 

The system Hamiltonian is split as Hs = Hq + H. 
The bath Hamiltonian is taken as Hb = 
where a\ and Ok are bosonic creation and annihilation 
operators respectively, ojk the frequency of each mode. 
The system-bath interaction is H[ = i(L (g) 61 — Ll 0 6), 
with b = where gk are coupling constants. The 

stationary bath density matrix ps is the vacuum state. 
In an interaction representation with respect to Hq and 
it follows 

Hi{t) = {iL)ZHt) + {-iL^)Z{t), (42) 

where the bath operator is Z{t) = '"^ith 

Qfe = ujk — uiQ- The frequency ujq follows from the assump¬ 
tion L{t) = e+^tHo ^ 



= -iH{t) Itpit)) + Zc{t)L \ip{t)) 



’ Sz,{s) ■ 


(43) 


From quantum measurement theory, this evolution can 
be read as a coherent unravelling of the bath (see Eq. 
(3.22) in [l^) where Z{t) is the noise operator. A hetero¬ 
dyne measurement process is recovered in the Markovian 
limit (Xlj . 

As an example of non-diagonal mapping, Eq. (I35L 
we take the single operator Z{t) = Z{t) -\- Z'^{t), and 
consistently with the diagonal case H {t) = Zq{t){iL). 
It is simple to check that Eq. (I38|) is fulfilled. In 
fact, given that Z^t) = Z{t), Eq. (|551) leads to the 
correlations Xq{t,s) = r]q{t,s) = TiB[pBZ{t)Z{s)] = 
Xfc (the same correlation than in the co¬ 

herent unraveling). Nevertheless, these correlations do 
not have associated a positive covariance matrix [Eq. 
(IA.4I) ]. This contradiction can be raised up by taking 
the complex noise Zq(t) as a real one, implying real cor¬ 
relations. Given that L ^ ^ here a real noise does 

not implies The previous complex correla¬ 

tions become real if we demand extra properties to the 
bath. If = —Qfe and, given that H\{t) — 
demanding g^k = gl, it follows Xq{i,s) = riq{t,s) = 
2Xfc>o l 5 fcPcos[Qfe(t - s)] = Pq{t - s). Therefore, from 
Eq. (fl^ we arrive to the generalized Schrodinger equa¬ 
tion 



= + Zq{t)L\tP{t)) 

-(Lt + L) f'dsPq{t-s)^-^^ 

Jo OZqlyS) 


(44) 


This evolution can be read as a q uadrature unravelling of 
the bath (see Eq. (4.30) in [ 23 ) where Z{t) Z^t) is 
the noise operator. A homodyne measurement process is 
recovered in the Markovian limit [ 2 ^ . 

We derived Eqs. (H51) and (1441) on the basis of the di¬ 
agonal and non-diagonal correlation mapping introduced 
previously. These results give a solid support to the 
present and previous analysis [ 2 ^ [ 2 ^, which rely on 
quantum measurement theory. 

Based on the analysis of Refs. [ 2 ^ [ 2 ^, we also con¬ 
clude that the ensemble of realizations associated to gen¬ 
eralized Schrodinger dynamics such as Eq. m cannot 
be read as conditional states of a system subjected to 









a continuous measurement process over the environment 
degrees of freedom. In fact, as argued in Refs. [H, [2^, 
in the non-Markovian regime the wave vector only repre¬ 
sents the state the system would be at a single time. The 
property of linking solutions at different times to make 
a trajectory of a continuously monitored system is lost 

M- 


V. SUMMARY AND CONCLUSIONS 

We based the present analysis on a stochastic ap¬ 
proach. Postulating an underlying evolution with 
multiplicative noises, we derived the generalized non- 
Markovian Gaussian stochastic Schrodinger equation 
(USD from the condition of average trace preservation. 
Complex Gaussian noises with arbitrary correlations 
s) and rjapit, s), Eqs. (|31) and (|1]) respectively, were 
considered. 

We focused our study on the symmetries of the ob¬ 
tained evolution. The generalized Schrodinger equation, 
after a redefinition of the noises, is invariant under arbi¬ 
trary unitary changes of the system operator base. This 
property allowed us to conclude that Eq. (HID and the 
stochastic dynamics derived in Ref. js^ are related by a 
unitary transformation. On the other hand, the associ¬ 
ated master equation share the same symmetry property 
than the stochastic Schrodinger equation. Therefore, in 
contrast to previous analysis with Markovian dynamics, 
we asked under which conditions the invariance property 
can be mapped with that of a microscopic (bosonian) 
system-bath dynamics, which in turn lead us to consider 
a mapping between the noise correlations Xa/ 3 it,s) and 
Vapity s) with the bath operator correlations. 

Two kind of maps were introduced. In the diagonal 
one, the quantum reservoir correlations are defined by the 
bath operators corresponding to the system-environment 
interaction Hamiltonian. In the non-diagonal map, a new 
set of bath operators is introduced (quadrature-like bath 
operators), being related to the previous ones by an ar¬ 
bitrary linear transformation. 

For the diagonal correlation mapping [Eq. (|5TD ]. the 
invariance symmetry property of the microscopic dynam¬ 
ics is shared by the stochastic unravelling only when the 
correlation r]afi{t,s) vanishes, Eqs. (1551) and (1551) . recov¬ 
ering in consequence the standard non-Markovian quan¬ 
tum diffusion model M- For the non-diagonal correla¬ 
tion map [Eq. (I36p ]. the invariance symmetry is shared 
only if some constraints are fulfilled. In fact, the inde¬ 
pendence of the density matrix evolution with respect to 
Va/sityS) and the positivity of the noise correlation ma¬ 
trix lead to the necessary conditions Eqs. (1551) and (IA.4I) 
respectively. These requirements establish the conditions 
under which the noise correlations that set the general¬ 
ized unraveling, Xajsitys) and riap{t,s), can be defined 
from the properties of the microscopic system-bath in¬ 
teraction. 

As example, we considered the case of a single noise 


channel defined by an Hermitian operator [Eq. (1391) ]. 
We concluded that only the standard version of the dy¬ 
namics can be related with a microscopic description if 
the invariance property is demanded. On the other hand, 
previous dynamics derived from quantum measurement 
theory were recovered from diagonal and non-diagonal 
correlation mapping [Eqs. (1551) and (|551) respectively] ap¬ 
plied to the same optical-like system-environment inter¬ 
action [2^ . These examples show that the consistence of 
the non-diagonal correlation map put severe constraints 
on the bath properties, which in general may or not be 
fulfilled. 

The present analysis not only clarifies the recent ad¬ 
vances in the formulation of Gaussian stochastic wave 
vector dynamics [3^ but also define the constraints under 
which the generalized unraveling can be put in one-to-one 
correspondence with a microscopic unitary description. 
They also supports the stringent constraints on the in¬ 
terpretation of non-Markovian Schrodinger equations as 
conditional measurement states at different times [^ . 
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Appendix: Arbitrary complex Ganssian noises 


The statistical properties of an arbitrary set of com¬ 
plex Gaussian noises {za{t)} can be defined through the 
characteristic functional 

G[/,5] = (|exp|i J dT[faiT)Za{T) + ga{T)z^{T)] 

(A.l) 

where / and g denotes a set of test functions {fa{t)} and 
{ga{t)}. The Gaussian statistics implies 



lnG[/,g] = 


- dr dsgaiT)xap{T,s)fp{s) 

Jo Jo 

-[ dr f dsfa{T)x*af)iT,s)gi3{s) 

Jo Jo 

- dr dsfa{T)r]al3{TyS)fi3{s) 

Jo Jo 

- f dr f dsgoy{T)glp{T,s)gp{s). (A.2) 
Jo Jo 


By functional derivation it follows that (zq, (t)) = 0, and 
the correlations defined by Eqs. ([SD and (HD. These defi¬ 
nitions cover the particular cases of real noises as well as 
the case where gap(t, s) = 0. 

Novikov’s theorem [371] gives an exact result for the 
mean value of the product between a Gaussian noise and 
any functional Ad of it. This theorem can be generalized 
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for the set of complex noises defined by Eq. (IA.2I) . Using 
the Gaussian property and the correlation definition it is 
possible to obtain 


{zy{t)M[{zc,{t)}])= dsx*p{t,s) 


w[{» W}] ' 

Szp{s) / 


+ 



dsri^p{t,s) 


5zi){s) r 


where M[{za{t)}] denotes a functional that depends on 
the set of noises {Za{t)}. The noise Zj{t) belongs to this 
set. This generalized Novikov theorem has been applied 
in the derivation of Sec. II. 


where n-yi^ is the average thermal number in each mode. 
Given that {t, s) = 0, the positivity constraint (IA.4I) 
(at any bath temperature) is satisfied whenever the ma¬ 
trixes of complex coefficients [{n^uj + l)ffa-y(w)g/ 37 (w)] 
and [n~fcuh’^^{uj)hp^{uj)] are positive defined, that is, 
the bath spectrum matrix is positive defined. On the 
other hand, we notice that non-stationary correlations 
[ 3 ^ Xap {t, s) ^ Xap {t ~ s), only arise if the under¬ 
lying Hamiltonian (EH) is time dependent. For exam¬ 
ple, the interaction strength may be time dependent, 
A Xip{t), Xap{t,s) ip^{t)xap jt - s ), or alterna¬ 
tively the previous interaction operators (IA.5I) are defined 
with time dependent coefficients gaa'{(^) 5 aa'(w,t), 

hacK'(^) ^ dctc^f {uJ^ . 


Correlation matrix 

The correlation matrix that define Eq. (|A.2I) has to 
satisfy a positivity constraint. By defining the scalar w = 
+ bi z'^{ti)), where a“ and 5“ are arbitrary 
coefficients, n € N, {U} arbitrary times, the absolute 
value > 0 implies that the four blocks correlation 
kernel 



must be positive defined, where x ^ XapiUjtj) and 
V VapiU^tj). An equivalent condition follows with 
w = + bflm[za{U)]), leading to 

( + + > 0 . (A.4b) 

V - v\ R-e[X -m J ~ 

In both cases, when n = I the positivity of the variance 
matrix of {za{t)} is recovered. Notice that r]ap(t,s) = 
Xa-pib, s) is an admissible correlation matrix only for real 
noises. 


a. Bosonic bath correlation matrix 

The correlation matrixes XaP {b, s) and rfap it, s) corre¬ 
sponding to the diagonal map, Eq. (ESI) and (ISil) respec¬ 
tively, satisfy the positivity constraint (IA.4I) . This prop¬ 
erty is fulfilled without imposing any special constraint 
on the bath properties or on the interaction Hamiltonian. 
In fact, by writing the interaction bath operators as lin¬ 
ear combinations of free bosonic modes 

Za{t) = f dw{gaa'iuj)b'l,^e'^^‘^* + haa'iuj)ba'u:e~"‘^*), 

(A.5) 

with commutation relations [bauj,b^p^,] = Sap6{u! — to'), it 
follows 

Xc.pit,s) = J duj{n^u> + 

+ J dujn^^hl^{uj)hpj{uj)e'^"‘^^*~'^\ (A.6) 


b. Exponential noise correlations 

As an (one dimensional) example, we assume that the 
realizations of the noise z{t) obey the linear stochastic 
differential equation 

^z{t) =-{x + in)z{t)+^{t), (A.7) 

where 7 > 0 and H are real parameters. The complex 

Gaussian white noise f (t) satisfies 

icmis)) = D6it-s), (A. 8 a) 

mm) = D'6{t-S), (A. 8 b) 

where D' < D are also real parameters. By integrating 
Eq. (IA.7I) as z{t) = z(0) exp[— ( 7 -|-zH)t]-|- f* dt' exp[— ( 7 -I- 
fH)(t — and assuming stationary initial condi¬ 

tions for z{t), it follows the exponential correlations 

X{t,s) = {z*{t)z{s)) = ^exp[-{x -in){t- s)], 

g{t,s) = {z{t)z{s)) = exp [-(7 -k iXl){t - s)], 

where t > s. In a Fourier domain, it is possible to demon¬ 
strate that these objects obey Eq. (IA.4I) . 

The degree of freedom introduced by ry(t, s) in Eq. (l39ll 
can be easily read by taking H = 0 , leading to 

X(t, s) - g{t, s) = (1 - ^ exp[- 7 (t - s)]. (A.9) 

Therefore, the dimensionless parameter D'/D allows to 
continuously departs from a Hamiltonian stochastic dy¬ 
namics {D' = D) and, in the other extreme, to reach the 
standard non-Markovian diffusion model, that is, D' = 0. 
While this property was derived for this particular case, 
it is simple to realize that the inclusion of the extra cor¬ 
relations r]ap(t,s) allow to smoothly reach these limits 
when considering the general evolution (1121) with an Her- 
mitian operator base, {L)^} = {La}. 
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